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We give an introduction to the heavy-quark effective theory and the 1/rriQ expan- 
sion, which provide the modern framework for a systematic, model-independent 
description of the properties and decays of hadrons containing a heavy quark. We 
discuss the applications of these concepts to spectroscopy and to the weak decays 
of B mesons. 



1 Introduction 

The weak decays of hadrons containing a heavy quark are employed for tests 
of the Standard Model and measurements of its parameters. They offer the 
most direct way to determine the weak mixing angles, to test the unitarity of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and to explore the physics 
of CP violation. At the same time, hadronic weak decays also serve as a probe 
of that part of strong-interaction phenomenology which is least understood: 
the confinement of quarks and gluons inside hadrons. 

The structure of weak interactions in the Standard Model is rather simple. 
Flavour-changing decays are mediated by the coupling of the charged current 
to the IF-boson field. At low energies, the charged-current interaction gives 
rise to local four-fermion couplings, whose strength is governed by the Fermi 
constant 

G F = / = 1.16639(2) GeV- 2 . (1) 

According to the structure of the these interactions, the weak decays of hadrons 
can be divided into three classes: leptonic decays, in which the quarks of the 
decaying hadron annihilate each other and only leptons appear in the final 
state; semileptonic decays, in which both leptons and hadrons appear in the 
final state; and non-leptonic decays, in which the final state consists of hadrons 
only. Representative examples of these three types of decays are shown in 
Fig. | 

The simple quark-line graphs shown in this figure are a gross oversim- 
plification, however. In the real world, quarks are confined inside hadrons, 
bound by the exchange of soft gluons. The simplicity of the weak interactions 
is overshadowed by the complexity of the strong interactions. A complicated 
interplay between the weak and strong forces characterizes the phenomenology 
of hadronic weak decays. As an example, a more realistic picture of a non- 
leptonic decay is shown in Fig. 0. Clearly, the complexity of strong-interaction 
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Figure 1: Examples of leptonic (B — > r u T ), semileptonic (B° — * D+e u e ), and non- 
leptonic (B° — > D+7r~) decays of B mesons. 



effects increases with the number of quarks appearing in the final state. Bound- 
state effects in leptonic decays can be lumped into a single parameter (a "de- 
cay constant"), while those in semileptonic decays are described by invariant 
form factors, depending on the momentum transfer q 2 between the hadrons. 
Approximate symmetries of the strong interactions help to constrain the prop- 
erties of these form factors. For non-leptonic decays, on the other hand, we are 
still far from having a quantitative understanding of strong-interaction effects 
even in the simplest decay modes. 

Over the last decade, a lot of information on heavy-quark decays has been 
collected in experiments at e + e~ and hadron colliders. This has led to a rather 
detailed knowledge of the flavour sector of the Standard Model and many of 
the parameters associated with it —There have been several great discoveries in 
this field, such as B°-B° mixing Era, charmless B decaysH ~0, and rare decays 
induced by penguin operators&Ll. The experimental progress in heavy-flavour 
physics has been accompanied by a significant progress in theory, which was 
related to the discovery of heavy-quark symmetry and the development of the 
heavy-quark effective theory (HQET). The excitement about these develop- 
ments is caused by the fact that they allow (some) model-independent predic- 
tions in an area in which "progress" in theory often meant nothing more than 
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Figure 2: More realistic representation of a non-leptonic decay. 



the construction of a new model, which could be used to estimate some strong- 
interaction hadronic matrix elements. Here we explain the physical picture 
behind heavy-quark symmetry and discuss the construction, as well as simple 
applications, of the heavy-quark expansion. Because of lack of time, we will 
have to focus on some particularly important aspects. A more complete dis- 
cussion of the applications of this formalismJjp heavy-flavour phenomenology 
can be found in some recent review artir-lpsPrM The reader is also encouraged 
to consult the earlier review papers EtlLj on the subject. 

Hadronic bound states of a heavy quark with light constituents (quarks, 
antiquarks and gluons) are characterized by a large separation of mass scales: 
the heavy-quark mass tuq is much larger than the mass scale Aqcd associated 
with the light degrees of freedom. Equivalently, the Compton wave length 
of the heavy quark (Aq ~ 1/rog) is much smaller than the size of the hadron 
containing the heavy quark (i?had ~ 1/Aqcd)- Our goal will be to separate the 
physics associated with these two scales, in such a way that all dependence on 
the heavy-quark mass becomes explicit. The framework in—which to perform 
this separation is the operator product expansion (OPE) LL3EZI. The HQET 
provides us with a convenient technical tool to construct the OPE. Before we 
start to explore the details of this effective theory, however, we should mention 
two important reasons why it is desirable to separate short- and long-distance 
physics in the first place: 

• A technical reason is that after the separation of short- and long-distance 
phenomena we can actually calculate a big portion of the relevant physics 
(i.e. all short-distance effects) using perturbation theory and renormali- 
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zation-group techniques. In particular, in this way we will be able to 
control all logarithmic dependence on the heavy-quark mass. 

• An important physical reason is that, after the short-distance physics has 
been separated, it may happen that the long-distance physics simplifies 
due to the realization of approximate symmetries, which imply non-trivial 
relations between observables. 

The second point is particularly exciting, since it allows us to make statements 
beyond the range of applicability of perturbation theory. Notice that here 
we are not talking about symmetries of the full QCD Lagrangian, such as its 
local gauge symmetry, but approximate symmetries realized in a particular 
kinematic situation. In particular, we will find that an approximate spin- 
flavour symmetry is realized in systems in which a single heavy quark interacts 
with light degrees of freedom by the exchange of soft gluons. 

At this point it is instructive to recall a more familiar example of how ap- 
proximate symmetries relate the long-distance physics of several observables. 
The strong interactions of pions are severely constrained by the approximate 
chiral symmetry of QCD. In a certain kinematic regime, where the momenta 
of the pions are much less than 1 GeV (the scale of chiral-symmetry breaking) , 
the long-distance physics of scattering amplitudes is encoded in a few "reduced 
matrix elements", such as the pion decay constant. An effective low-energy 
theory called chiral perturbation theory provides a systematic expansion of 
scattering amplitudes in powers of the pion momenta, and thus helps to de- 
rive the relations between different scattering amplitudes imposed by chiral 
symmetry lia. We will find that a similar situation holds for the case of heavy 
quarks. Heavy-quark symmetry implies that, in the limit where itiq 3> Aqcd, 
the long-distance physics of several observables is encoded in few hadronic pa- 
rameters, which can be defined in terms of operator matrix elements in the 
HQET. 

2 Heavy- Quark Symmetry 

2.1 The Physical Picture 

There are several reasons why the strong interactions of systems contain- 
ing heavy quarks are easier to understand than those of systems containing 
only light quarks. The first is asymptotic freedom, the fact that the effective 
coupling constant of QCD becomes weak in processes with la rge m omentum 
transfer, corresponding to interactions at short-distance scales Ej'El At large 
distances, on the other hand, the coupling becomes strongcl, leading to non- 
perturbative phenomena such as the confinement of quarks and gluons on a 



4 



length scale i?had ~ 1/Aqcd ~ 1 &n ; which determines the size of hadrons. 
Roughly speaking, Aqcd ~ 0.2 GeV is the energy scale that separates the 
regions of large and small coupling constant. When the mass of a quark Q 
is much larger than this scale, it is called a heavy quark. The quarks of the 
Standard Model fall naturally into two classes: up, down and strange are light 
quarks, whereas charm, bottom and top are heavy quarks.^] For heavy quarks, 
the effective coupling constant a s (mQ) is small, implying that on length scales 
comparable to the Compton wavelength Aq ~ 1/mg the strong interactions 
are perturbative and similar to the electromagnetic interactions. In fact, the 
quarkonium systems {QQ), whose size is of order Aq/a s (mq) <C i?had, are 
very much hydrogen-like. 

Systems composed of a heavy quark and light constituents are more com- 
plicated, however. The size of such systems is determined by i?had, and the 
typical momenta exchanged between the heavy and light constituents are of 
order Aqcd- The heavy quark is surrounded by a most complicated, strongly 
interacting cloud of light quarks, antiquarks, and gluons. In this case it is the 
fact that Aq <C -Rhad, i.e. that the Compton wavelength of the heavy quark is 
much smaller than the size of the hadron, which leads to simplifications. To 
resolve the quantum numbers of the heavy quark would require a hard probe; 
the soft gluons exchanged between the heavy quark and the light constituents 
can only resolve distances much larger than Aq. Therefore, the light degrees 
of freedom are blind to the flavour (mass) and spin orientation of the heavy 
quark. They experience only its colour field, which extends over large dis- 
tances because of confinement. In the rest frame of the heavy quark, it is in 
fact only the electric colour field that is important; relativistic effects such as 
colour magnetism vanish as mq — > oo. Since the heavy-quark spin partici- 
pates in interactions only through such relativistic effects, it decouples. That 
the heavy-quark mass becomes irrelevant can be seen as follows: As mq — » oo, 
the heavy quark and the hadron that contains it have the same velocity. In the 
rest frame of the hadron, the heavy quark is at rest, too. The wave function 
of the light constituents follows from a solution of the field equations of QCD 
subject to the boundary condition of a static triplet source of colour at the 
location of the heavy quark. This boundary condition is independent of mq, 
and so is the solution for the configuration of the light constituents. 

It follows that, in the limit mq — > oo, hadronic systems which differ only 
in the flavour or spin quantum numbers of tha-heavy quark have the same 
configuration of their light degrees of freedom E3~E3. Although this observa- 
tion still does not allow us to calculate what this configuration is, it provides 

"Ironically, the top quark is of no relevance to our discussion here, since it is too heavy to 
form hadronic bound states before it decays. 



5 



relations between the properties of such particles as the heavy mesons B, D, 
B* and D* , or the heavy baryons A& and A c (to the extent that corrections 
to the infinite quark-mass limit are small in these systems). These relations 
result from some approximate symmetries of the effective strong interactions 
of heavy quarks at low energies. The configuration of light degrees of freedom 
in a hadron containing a single heavy quark with velocity v does not change 
if this quark is replaced by another heavy quark with different flavour or spin, 
but with the same velocity. Both heavy quarks lead to the same static colour 
field. For Nh heavy-quark flavours, there is thus an SU^A^) spin-flavour sym- 
metry group, under which the effective strong interactions are invariant. These 
symmetries are in close correspondence to familiar properties of atoms: The 
flavour symmetry is analogous to the fact that different isotopes have the same 
chemistry, since to a good approximation the wave function of the electrons 
is independent of the mass of the nucleus. The electrons only see the total 
nuclear charge. The spin symmetry is analogous to the fact that the hyperfine 
levels in atoms are nearly degenerate. The nuclear spin decouples in the limit 
m e /mN — > 0. 

Heavy-quark symmetry is an approximate symmetry, and corrections arise 
since the quark masses are not infinitely heavy. In many respects, it is comple- 
mentary to chiral symmetry, which arises in the opposite limit of small quark 
masses. However, whereas chiral symmetry is a symmetry of the QCD La- 
grangian in the limit of vanishing quark masses, heavy-quark symmetry is not 
a symmetry of the Lagrangian (not even an approximate one), but rather a 
symmetry of an effective theory, which is a good approximation of QCD in a 
certain kinematic region. It is realized only in systems in which a heavy quark 
interacts predominantly by the exchange of soft gluons. In such systems the 
heavy quark is almost on shell; its momentum fluctuates around the mass shell 
by an amount of order Aqcd- The corresponding fluctuations in the velocity of 
the heavy quark vanish as Aqcd/toq — > 0. The velocity becomes a conserved 
quantity and is no longer a dynamical degree of freedom E3. Nevertheless, re- 
sults derived on the basis of heavy-quark symmetry are model-independent 
consequences of QCD in a well-defined limit. The symmetry-breaking correc- 
tions can, at least in principle, be studied in a systematic way. A convenient 
framework for analysing these corrections is provided by the HQET. Before 
presenting a detailed discussion of the formalism of this effective theory, we 
shall first point out some of the important implications of heavy-quark sym- 
metry for the spectroscopy and weak decays of heavy hadrons. 
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2.2 Spectroscopic Implications 

The spin-flavour symmetry leads to many interesting relations between the 
properties of hadrons containing a heavy nuark. The most direct consequences 
concern the spectroscopy of such states E2I. In the limit mg — ► oo, the spin 
of the heavy quark and the total angular momentum j of the light degrees of 
freedom inside a hadron are separately conserved by the strong interactions. 
Because of heavy-quark symmetry, the dynamics is independent of the spin and 
mass of the heavy quark. Hadronic states can thus be classified by the quantum 
numbers (flavour, spin, parity, etc.) of the light degrees of freedom El. The 
spin symmetry predicts that, for fixed j ^ 0, there is a doublet of degenerate 
states with total spin J = j ± |. The flavour symmetry relates the properties 
of states with different heavy-quark flavour. 

In general, the mass of a hadron Hq containing a heavy quark Q obeys an 
expansion of the form 

Am 2 

m„=m Q +A+- +0(l/m%). (2) 

2m Q w 

The parameter A represents contributions arising frorn_all terms in the La- 
grangian that are independent of the heavy-quark massu, whereas the quan- 
tity Am 2 originates from the terms of order 1 / ttlq in the effective Lagrangian 
of the HQET. For the moment, the detailed structure of these terms is of no 
relevance; it will be discussed at length later on. For the ground-state pseu- 
doscalar and vector mesons, one can parametrize the contributions from the 
1/toq corrections in terms of two quantities, Ai and A2, in such a way thattj 



Am 2 = -Ai + 2 



J(J + 1)-| A 2 . (3) 



Here J is the total spin of the meson. The first term, — Ai/2mQ, arises from the 
kinetic energy of the heavy quark inside the meson; the second term describes 
the interaction of the heavy-quark spin with the gluon field. The hadronic pa- 
rameters A, Ai and A2 are independent of toq. They characterize the properties 
of the light constituents. 

Consider, as a first example, the SU(3) mass splittings for heavy mesons. 
The heavy-quark expansion predicts that 

m Bs - m Bd = A s - A d + 0(l/m b ) , 

m Ds - m Dd = A s - A d + 0(1 /m c ) , (4) 

where we have indicated that the value of the parameter A depends on the 
flavour of the light quark. Thus, to the extent that the charm and bottom 
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quarks can both be considered sufficiently heavy, the mass splittings should 
be similar in the two systems. This prediction is confirmed experimentally, 
since c3 

m Bs -m Bd = (90 ±3) MeV, 

m Ds - m Dd = (99 ± 1) MeV . (5) 

As a second example, consider the spin splittings between the ground-state 
pseudoscalar (J = 0) and vector (J = 1) mesons, which are the members of 
the spin-doublet with j = |. The theory predicts that 

m%, - m% = 4A 2 + 0(l/m b ) , 

m 2 D , - m 2 D = 4A 2 + 0(l/m c ) . (6) 
The data are compatible with this: 

m|. - m% ~ 0.49 GeV 2 , 

m 2 D * -m 2 D ~ 0.55 GeV 2 . (7) 

Assuming that the B system is close to the heavy-quark limit, we obtain the 
value 

A 2 ~ 0.12 GeV 2 (8) 

for one of the hadronic parameters in (^|). As we shall see later, this quantity 
plays an important role in the phenomenology of inclusive decays of heavy 
hadrons. 

A third example is provided by the mass splittings between the ground- 
state mesons and baryons containing a heavy quark. The HQET predicts that 

m\ b - m B = Abaryon - A mcson + 0(l/m b ) , 

m Ac - m D = Abaryon - A mcson + 0(1/to c ) . (9) 

This is again consistent with the experimental results 

toa„ - m B = (346 ± 6) MeV , 

m Ac - m D = (416 ± 1) MeV , (10) 

although in this case the data indicate sizeable symmetry-breaking corrections. 
For the mass of the At baryon, we have used the value 

m A b = (5625 ± 6) MeV , (11) 
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which is obtained by averaging the result E3 m Ab = (5639 ± 15) MeV with 
the value m Ab = (5623 ± 5 ± 4) MeV reported by the CDF Collaboration E3. 
The dominant correction to the relations ([)]) comes from the contribution of 
the chromo-magnetic interaction to the masses of the heavy mesons,^ which 
adds a term ZXifluiQ on the right-hand side. Including this term, we obtain 
the refined prediction that the values of the following two quantities should be 
close to each other: 

m Ab - m B - — 1 = (312 ± 6) MeV , 
2m B 

o \ 

m\ c -m D = (320 ± 1) MeV (12) 

2m rj 

This is clearly satisfied by the data. 

The mass formula (^|) can also be used to derive information on the heavy- 
quark (pole) masses from the observed hadron masses. Introducing the "spin- 
averaged" meson masses tob = \ (ms+3ms*) ~ 5.31 GeVandmu = \ (mo + 
3mo«) — 1.97 GeV, we find that 

m h - m c = (M B - m D ) { 1 - n _^ 1 _ + 0(l/m%) \ , (13) 
[ 2m B m D ^ J 

where 0(1/itiq) is used as a generic notation representing terms suppressed 
by three powers of the b- or c-quark masses-. Using theoretical estimates for 
the parameter Ai, which lie in the rangeE3 l3 

Ai = -(0.3 ±0.2) GeV 2 , (14) 

this relation leads to 



77jf, — m c — (3.39 i 0.03 ± 0.03) GeV , (15) 

where the first error reflects the uncertainty in the value of Ai, and the second 
one takes into account unknown higher-order corrections. 



2.3 Exclusive Semileptonic Decays 

Semileptonic decays of B mesons have received a lot of attention in recent 
years. The decay channel B — > D*£ v has the largest branching fraction of all 
-B-meson decay modes. From a theoretical point of view, semileptonic decays 
are simple enough to allow for a reliable, quantitative description. The analysis 

^Because of the spin symmetry, there is no such contribution to the masses of the Aq baryons. 
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of these decays provides much information about the strong forces that bind 
the quarks and gluons into hadrons. Heavy-quark symmetry implies relations 
between the weak decay form factors of heavy mesons, which are of particular 
interest. These relations have been derived by Isgur and WiseEH, generalizing 
ideas developed by Nussinov and WetzelE.3, and by Voloshin and ShifmanE3tJ. 

Consider the elastic scattering of a B meson, B(v) — ► -B(u'), induced by 
a vector current coupled to the b quark. Before the action of the current, the 
light degrees of freedom inside the B meson orbit around the heavy quark, 
which acts as a static source of colour. On average, the b quark and the 
B meson have the same velocity v. The action of the current is to replace 
instantaneously (at time t = to) the colour source by one moving at a velocity 
v' , as indicated in Fig. ||. If v = v', nothing happens; the light degrees of 
freedom do not realize that there was a current acting on the heavy quark. 
If the velocities are different, however, the light constituents suddenly find 
themselves interacting with a moving colour source. Soft gluons have to be 
exchanged to rearrange them so as to form a B meson moving at velocity v' . 
This rearrangement leads to a form-factor suppression, which reflects the fact 
that as the velocities become more and more different, the probability for an 
elastic transition decreases. The important observation is that, in the limit 
rat — > co, the form factor can only depend on the Lorentz boost 7 = v • v' 
that connects the rest frames of the initial- and final-state mesons. Thus, in 
this limit a dimensionless probability function £(u • v') describes the transition. 
It is called the Isgur-Wise function E3. In the HQET, which provides the 
appropriate framework for taking the limit m& — * 00, the hadronic matrix 
element describing the scattering process can thus be written as 

— (B(v')\ bv^by \B(v)) = £(v ■ v') {v + v'Y ■ (16) 

Here, b v and b v > are the velocity-dependent heavy-quark fields of the HQET, 
whose precise definition will be discussed later. It is important that the func- 
tion £(d ■ v') does not depend on mi,. The factor l/mg on the left-hand side 
compensates for a trivial dependence on the heavy-meson mass caused by the 
rclativistic normalization of meson states, which is conventionally taken to be 

(B(p')\B(p)) = 2m B v (2nfS 3 (p-p'). (17) 

Note that there is no term proportional to (v — v')^ in ([l6[). This can be seen 
by contracting the matrix element with (v — v')^, which must give zero since 
j)b v = b v and b V 'ft = b v >. 

It is more conventional to write the above matrix element in terms of an 
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Figure 3: Elastic transition induced by an external heavy-quark current. 

clastic form factor F c \(q 2 ) depending on the momentum transfer q 2 = (p—p') 2 : 

(B(v')\h%\B(v)) = F cl (q 2 )(p + P y , (18) 

where p^ = mBV^'\ Comparing this with (|l6|), we find that 

F cl (q 2 ) = £(v ■ v') , q 2 = -2m%(v ■ v' - 1) . (19) 

Because of current conservation, the elastic form factor is normalized to unity 
at q 2 = 0. This condition implies the normalization of the Isgur-Wise function 
at the kinematic point v ■ v' = 1, i.e. for v = v': 

£(!)=!■ (20) 

It is in accordance with the intuitive argument that the probability for an 
clastic transition is unity if there is no velocity change. Since for v — v' the 
daughter meson is at rest in the rest frame of the parent meson, the point 
v ■ v' = 1 is referred to as the zero-recoil limit. 

We can now use the flavour symmetry to replace the b quark in the final- 
state meson by a c quark, thereby turning the B meson into a D meson. Then 
the scattering process turns into a weak decay process. In the infinite-mass 
limit, the replacement b v > — > c v > is a symmetry transformation, under which 
the effective Lagrangian is invariant. Hence, the matrix element 

' 'D(v')\ c v ,rb v \B(v)) = £(v ■ v') (v + v'Y (21) 



is still determined by the same function £(v ■ v'). This is interesting, since in 
general the matrix element of a flavour-changing current between two pseu- 
doscalar mesons is described by two form factors: 

(D(v')\ c^b \B{v)) = f+{q 2 ) (p + p'Y - f-(<l 2 ) (P - P'Y ■ (22) 
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Comparing the above two equations, we find that 



2 m B ± m D 
2^/m B m D 

q 2 = m B + m 2 D — 2vn B m B> v ■ v' . (23) 

Thus, the heavy-quark flavour symmetry relates two a priori independent form 
factors to one and the same function. Moreover, the normalization of the 
Isgur-Wise function at v ■ v' = 1 now implies a non-trivial normalization of 
the form factors f±(q 2 ) at the point of maximum momentum transfer, g^ ax = 
(m B - m D ) 2 : 

2 m B ±m D 

/±(gmax) = 9 - = ■ ( 24 ) 

The heavy-quark spin symmetry leads to additional relations among weak 
decay form factors. It can be used to relate matrix elements involving vector 
mesons to those involving pseudoscalar mesons. A vector meson with longi- 
tudinal polarization is related to a pseudoscalar meson by a rotation of the 
heavy-quark spin. Hence, the spin-symmetry transformation cjj, — > c^-xelates 
B — ► D with B — ► D* transitions. The result of this transformation is Ell: 

' (D*(v',e)\c v ,^b v \B(v)) = ie^s*^ £(«■«'), 



y<m B mD* 

! (D*{v',e)\c v ^ 5 K\B(v)) = s** (v ■ v' + 1) - v'" e* ■ v £(v ■ v') 



y/m B m D * 



(25) 



where e denotes the polarization vector of the D* meson. Once again, the 
matrix elements are completely described in terms of the Isgur-Wise function. 
Now this is even more remarkable, since in general four form factors, V(q 2 ) for 
the vector current, and Ai(q 2 ), i — 0,1,2, for the axial current, are required 
to parametrize these matrix elements. In the heavy-quark limit, they obey the 
relations 



- V(q 2 )=A (q 2 )=A 1 (q 2 ) 



1 q 



2 



-1 



(m B + m D ) 2 
j 2 = m 2 B + m 2 D , — 2m B mn* v ■ v' . (26) 
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Equations (|23|) and ( |26| ) summarize the relations imposed by heavy-quark 
symmetry on the weak decay form factors describing the semileptonic decay 
processes B — > D Iv and B — > D*lv. These relations are model-independent 
consequences of QCD in the limit where mi,, m c Aqcd- They play a crucial 
role in the determination of the CKM matrix element | V c b | ■ In terms of the 
recoil variable w — V-^>\ the differential semileptonic decay rates in the heavy- 
quark limit becomet3: 

6r{B 7w Di9) = i& lVcb]2 {mB + mD? ml {w2 ~ 1)3/2 e[w) ' 



dw 487T 3 
1 + 



iVcbl^ms-TOc-) 2 ™!,, a/w 2 - 1 (iu + 1) 

Aw m B — 2wmsmD* + fn^i 
w + 1 (tub — rriD*) 2 



rM-(27) 



These expressions receive symmetry-breaking corrections, since the masses of 
the heavy quarks are not infinitely heavy. Perturbative corrections of order 
a s( m Q) can be calculated order by order in perturbation theory. A more 
difficult task is to control the non-perturbative power corrections of order 
(Aqcd/wq)™- The HQET provides a systematic framework for analysing these 
corrections. For the case of weak-decay form factors, the analysis of the 1/toq 
corrections was performed by LukeCj. Later, Falk and the present author have 
also analysed the structure of 1/tjIq corrections for both meson and baryon 
weak decay form factors t3. We shall not discuss these rather technical issues 
in detail, but only mention the most important result of Luke's analysis. It 
concerns the zero- recoil limit, where an analogue of the Ademollo-Gatto theo- 
remlid can be proved. This is Luke's theoremLJ, which states that the matrix 
elements describing the leading X/mq corrections to weak decay amplitudes 
vanish at zero recoil. This theorem is valid to all orders in perturbation the- 
ory EjujH Most importantly, it protects the B —±-I)*£D decay rate from 
receiving first-order 1/toq corrections at zero recoil E£l. (A similar statement 
is not true for the decay B — ► D £ 9, however. The reason is simple but some- 
what subtle. Luke's theorem protects only those form factors not multiplied by 
kinematic factors that vanish for v = v' . By angular momentum conservation, 
the two pseudoscalar mesons in the decay B — > D iv must be in a relative p 
wave, and hence the amplitude is proportional to the velocity I^dI of the D 
meson in the B-meson rest frame. This leads to a factor (w 2 — 1) in the decay 
rate. In such a situation, form factors that are kinematically suppressed can 
contributeE3.) 
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2.4 Model-Independent Determination of \V c b\ 



We will now discuss the most important application of the HQET in the context 
of semileptonic decays of B mesons. A model-independent determination of the 
CKM matrix element | V c b \ based on heavy-quark symmetry can be obtainedJay 
measuring the recoil spectrum of D* mesons produced in B — ► D*£ v decaysEj. 
In the heavy-quark limit, the differential decay rate for this process has been 
given in (^7|). In order to allow for corrections to that limit, we write 

dT(B^D*£v) _ G 



dw 



'' (m B — m D *) 2 m 3 D , \/ w 2 — 1 (w + 1) 
2w msmD* + fn 



1 + 



Aw m 2 B 



D> 



w 



1 



(m B - m D * 



\V cb \ 2 T 2 (w), 
(28) 



where the hadronic form factor T(w) coincides with the Isgur-Wise function 
up to symmetry-breaking corrections of order a s (m,Q) and Aqcd/wq. The 
idea is to measure the product \V c b\ 3~(w) as a function of w, and to extract 
\Vcb\ from an extrapolation of the data to the zero-recoil point w = 1, where 
the B and the D* mesons have a common rest frame. At this kinematic point, 
heavy-quark symmetry helps to calculate the normalization with small 

and controlled theoretical errors. Since the range of w values accessible in this 
decay is rather small (1 < w < 1.5), the extrapolation can be done using an 
expansion around w = 1: 



T{w) = 1 - Q 2 (w - 1 



(29) 



The slope g 2 is treated as a fit parameter. 

Measurements of_the recoil spectrum have been performed first by the 
ARGUSo and CLEOllj Collaborations in experiments operatixig at the T(4s) 
resonance, and more recently by the ALEPHC3 and DELPHIcZI Collaborations 
at LEP. As an example, Fig. || shows the data reported by the CLEO Collabo- 
ration. The results obtained by the various experimental groups from a linear 
fit to their data are summarized in Table |]. The weighted average of these 
results is 



\V cb \T(l) = (34.6 ±1.7) 
g 2 = 0.82 ±0.09. 



x 10 



-3 



(30) 



The ef 
Stone c 



2ct of a positive curvature of the form factor has been investigated by 
J, who finds that the value of \V c b\ ^"(1) may change by up to +4%. We 
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Figure 4: CLEO rtaita for the product \V c b\ F(w), as extracted from the recoil spectrum in 
B — > D*lv decaystij. The line shows a linear fit to the data. 

thus increase the above value by (2 ± 2)% and quote the final result as 

\V ch \T{\) = (35.3 ± 1.8) x 1CT 3 . (31) 

In future analyses, the extrapolation to zero recoil should be performed in- 
cluding higher-order terms in the expansion (^9|) . It can be shown in a model- 
independent way that the shape of the fnrrn factor is highly constrained by 
analyticity and unitarity requirements Ej'lj. In particular, the curvature at 
w = 1 is strongly correlated with the slope of the form factor. For the value of 
g 2 given in (|30|), one obtains a small positive curvature Ej, in agreement with 
the assumption made in Ref. 48. 

Table 1: Values for V^j, | ^-"(1) (in units of 10 — 3 ) and § 2 extracted from measurements of 
the recoil spectrum in B — > D*lv decays 





\v ch \m (10- 3 ) 


Q 2 


ARGUS 
CLEO 
ALEPH 
DELPHI 


38.8 ±4.3 ±2.5 
35.1 ±1.9 ±2.0 
31.4 ±2.3 ±2.5 
35.0 ± 1.9 ±2.3 


1.17 ±0.22 ±0.06 
0.84 ±0.12 ±0.08 
0.39 ±0.21 ±0.12 
0.81 ±0.16 ±0.10 



Heavy-quark symmetry implies that the general structure of the symmetry- 
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breaking corrections to the form factor at zero recoil isLJ 

jr(l) = 7 1A ( 1 + x + const x ^QCD + \ ^ (1 + tf j (32) 

V m Q m Q J 

where 7?a is a short-distance correction arising from the (finite) renormalization 
of the flavour-changing axial current at zero recoil, and 5\/ m 2 parametrizes 
second-order (and higher) power corrections. The absence of -first-order power 
corrections at zero recoil is a consequence of Luke^s-.thenremcj. The one-loop 
expression for t/a has been known for a long timeE3i23'Ell: 

+ as(M) ( rnt + rnc ^ n* _ 8\ ^ q gg (33) 
7r \ nit — M c m c 3 / 

The scale M in the running coupling constant can bejixed by adopting the pre- 
scription of Brodsky, Lepage and Mackenzie (BLM)E3, according to which it is 
identified with the average virtuality of the gluon in the one- loop diagrams that 
contribute to if a- If QulM) is defined in the modified minimal subtraction (Ms) 
scheme, the result isE3 M ~ 0.51^/r7i c mb. Several estimates of higher-order 
corrections to r]A have been discussed. The next-to-l eadin g order resummation 
of logarithms of the type [a s ln(mb/m c )]" leads tooE3 t\a — 0.985. On the 
other hand, the resummation of "renormalon-chain" contributions of the form 
/3o _1 a™, where /3q is the first coefficient of the QCD /3-function, gives l£k3 
r\A — 0.945. Using these partial resummations to estimate the uncertainty 
results in t]a = 0.965 ± 0.020. Recently, Czarnecki has improved this estimate 
by calculating t]a at two-loop orderES. His result, 

r)A = 0.960 ± 0.007 , (34) 

is in excellent agreement with the BLM- improved one-loop estimate (^). Here 
the error is taken to be the size of the two-loop correction. 

The analysis of the power corrections 8i/ m 2 is more difficult, since it cannot 
rely on perturbation theory. Three approaches have been discussed: in the 
"exclusive approach" , all 1/ nig operatprsin the HQET are classified and their 
matrix elements estimated, leading toc3H Si/ m 2 — —(3 ± 2)%; the "inclusive 
approach" has been used to derive the bound 5i/ m 2 < —3%, and to estimate 
thatEr[] 6x/ m 2 — —(7 ± 3)%; the "hybrid approach" combines the virtues of 
thc_formcr two to obtain a more restrictive lower bound on <$i/ m 2. This leads 

toH 

Si/m* = -0.055 ±0.025. (35) 

c This bound has been criticised in Ref. 61. 
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Combining the above results, adding the theoretical errors linearly to be 
conservative, gives 

= 0.91 ±0.03 (36) 

for the normalization of the hadronic form factor at zero recoil. Thus, the 
corrections to the heavy-quark limit amount to a moderate decrease of the 
form factor of about 10%. This can be used to extract from the experimental 
result ( [31] ) the model-independent value 

\Vcb\ = (38.8 ± 2.0 cxp ± 1.2th) x 10~ 3 . (37) 

3 Heavy- Quark Effective Theory 

3.1 The Effective Lagrangian 

The effects of a very heavy particle often become irrelevant at low energies. It 
is then useful to construct pa, low-energy effective theory, in which this heavy 
particle no longer appears E3. Eventually, this effective theory will be easier 
to deal with than the full theory. A familiar example is Fermi's theory of 
the weak interactions. For the description of weak decays of hadrons, the 
weak interactions can be approximated by point- like four-fermion couplings, 
governed by a dimcnsionful coupling constant Gf- Only at energies much 
larger than the masses of hadrons can the effects of the intermediate vector 
bosons, W and Z, be resolved. 

The process of removing _the degrees of freedom of a heavy particle in- 
volves the following steps E-aE3: one first identifies the heavy-particle fields 
and "integrates them out" in the generating functional of the Green functions 
of the theory. This is possible since at low energies the heavy particle does not 
appear as an external state. However, although the action of the full theory 
is usually a local one, what results after this first step is a non-local effective 
action. The non-locality is related to the fact that in the full theory the heavy 
particle with mass M can appear in virtual processes and propagate over a 
short but finite distance Ax ~ 1/M. Thus, a second step is required to obtain 
a local effective Lagrangian: the non-local effective action is rewritten as an 
infinite series of local terms using the OPEEjEJ. Roughly speaking, this corre- 
sponds to an expansion in powers of 1/M. It is in this step that the short- and 
long-distance physics is disentangled. The long-distance physics corresponds 
to interactions at low energies and is the same in the full and the effective 
theory. But short-distance effects arising from quantum corrections involving 
large virtual momenta (of order M) are not reproduced in the effective theory 
once the heavy particle has been integrated out. In a third step, they have to 
be added in a perturbative way using renormalization-group techniques. These 
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short-distance effects lead to a renormalization of the coefficients of the local 
operators in the effective Lagrangian. An example is the effective Lagrangian 
for non-leptonic weak decays, in which radiative corrections from hard gluons 
with virtual momenta in the range between M^v and some renormalization 
scale /j, ~ 1 GeV give ris£_.to_Wilson coefficients, which renormalize the local 
four-fermion interactionsEJlia. 
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Figure 5: Philosophy of the heavy-quark effective theory. 



The heavy-quark effective theory (HQET) is constructed to provide a sim- 
plified description of processes where a heavy quark interact^with light degrees 
of freedom predominantly by the exchange of soft gluons L2l Clearly, itlq 
is the high-energy scale in this case, and Aqcd is the scale of the hadronic 
physics we are interested in. The situation is illustrated in Fig. ^. At short 
distances, i.e. for energy scales larger than the heavy-quark mass, the physics is 
perturbative and described by ordinary QCD. For mass scales much below the 
heavy-quark mass, the physics is complicated and non-perturbative because 
of confinement. Our goal is to obtain a simplified description in this region 
using an effective field theory. To separate short- and long-distance effects, we 
introduce a separation scale \i such that Aqcd </J< vtiQ . The HQET will be 
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constructed in such a way that it is identical to QCD in the long-distance re- 
gion, i.e. for scales below /i. In the short-distance region, the effective theory is 
incomplete, however, since some high-momentum modes have been integrated 
out from the full theory. The fact that the physics must be independent of the 
arbitrary scale ji allows us to derive renormalization-group equations, which 
can be employed to deal with the short-distance effects in an efficient way. 

Compared with most effective theories, in which the degrees of freedom 
of a heavy particle are removed completely from the low-energy theory, the 
HQET is special in that its purpose is to describe the properties and decays of 
hadrons which do contain a heavy quark. Hence, it is not possible to remove 
the heavy quark completely from the effective theory. What is possible is to 
integrate out the "small components" in the full heavy-quark spinor, which 
describe the fluctuations around the mass shell. 

The starting point in the construction of the low-energy effective theory 
is the observation that a very heavy quark bound inside a hadron moves more 
or less with the hadron's velocity v, and is almost on shell. Its momentum can 
be written as 

p£ = m Q v» + , (38) 

where the components of the so-called residual momentum k are much smaller 
than toq . Note that v is a four- velocity, so that v 2 = 1 . Interactions of the 
heavy quark with light degrees of freedom change the residual momentum by 
an amount of order Afc <~ Aqcd, but the corresponding changes in the heavy- 
quark velocity vanish as Aqcd/tiq — ► 0. In this situation, it is appropriate to 
introduce large- and small-component fields, h v and H v , by 

h v (x) = e lm ^ x P+ Q{x) , H v (x) = e im ^ x P_ Q(x) , (39) 

where P+ and P_ are projection operators defined as 

P± = ±¥- (40) 

It follows that 

Q(x) = e - im * v - x [h v (x) + H v (x)} . (41) 

Because of the projection operators, the new fields satisfy f h v = h v and f H v = 
—H v . In the rest frame, i.e. for = (1,0,0,0), h v corresponds to the upper 
two components of Q, while H v corresponds to the lower ones. Whereas h v 
annihilates a heavy quark with velocity v, H v creates a heavy antiquark with 
velocity v. 
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In terms of the new fields, the QCD Lagrangian for a heavy quark takes 
the form 

£q = Q(ip- tuq) Q 

— h v iv-D h v — H v (iv ■ D + 2rriQ ) H v 

+ h v ip^H v +H v ip^h v , (42) 

where Dj_ = — v^v-Dis orthogonal to the heavy-quark velocity: v ■ D± = 0. 
In the rest frame, Dj = (0, D ) contains the spatial components of the covariant 
derivative. From ((42|), it is apparent that h v describes massless degrees of 
freedom, whereas H v corresponds to fluctuations with twice the heavy-quark 
mass. These are the heavy degrees of freedom that will be eliminated in the 
construction of the effective theory. The fields are mixed by the presence of the 
third and fourth terms, which describe pair creation or annihilation of heavy 
quarks and antiquarks. As shown in the first diagram in Fig. |, in a virtual 
process a heavy quark propagating forward in time can turn into an antiquark 
propagating backward in time, and then turn back into a quark. The energy 
of the intermediate quantum state hhH is larger than the energy of the initial 
heavy quark by at least 2m,Q. Because of this large energy gap, the virtual 
quantum fluctuation can only propagate over a short distance Ax ~ 1/ttiq. 
On hadronic scales set by i?had ~ 1/Aqcd, the process essentially looks like a 
local interaction of the form 

h v ip ± -^—ip x h v , (43) 
2m Q 

where we have simply replaced the propagator for H v by l/2m<g. A more 
correct treatment is to integrate out the small-component field H v , thereby 
deriving a non-local effective action for the large-component field h v , which 
can then be expanded in terms of local operators. Before doing this, let us 
mention a second type of virtual corrections involving pair creation, namely 
heavy-quark loops. An example is shown in the second diagram in Fig. ^[ 
Heavy-quark loops cannot be described in terms of the effective fields h v and 
H v , since the quark velocities inside a loop are not conserved and are in no 
way related to hadron velocities. However, such short-distance processes are 
proportional to the small coupling constant a s (raq) and can be calculated in 
perturbation theory. They lead to corrections that are added onto the low- 
energy effective theory in the renormalization procedure. 

On a classical level, the heavy degrees of freedom represented by H v can 
be eliminated using the equation of motion. Taking the variation of the La- 
grangian with respect to the field H Vl we obtain 

(ivD + 2m Q )H v = ip ± h v . (44) 
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Figure 6: Virtual fluctuations involving pair creation of heavy quarks. In the first diagram, 
time flows to the right. 



This equation can formally be solved to give 

2tiiq + iv-D 

showing that the small-component field H v is indeed of order 1/wq. We can 
now insert this solution into (fl2|) to obtain the "non-local effective Lagrangian" 

C cS = h v iv-D h v + h v i p± — — —ip±h v . (46) 

IrtiQ + iv- D 

Clearly, the second term corresponds to the first class of virtual processes 
shown in Fig. ^. 

It is possible to derive this Lagrangian in a more elegant way by manipu- 
lating the generating functional for QCD Green's functions containing heavy- 
quark fields E3. To this end, one starts from the field redefinition ( |4l| ) and 
couples the large-component fields h v to external sources p v . Green's functions 
with an arbitrary number of h v fields can be constructed by taking derivatives 
with respect to p v . No sources are needed for the heavy degrees of freedom 
represented by H v . The functional integral over these fields is Gaussian and 
can be performed explicitly, leading to the effective action 



S* cff = / d 4 x£ cff -ilnA, (47) 

with C e s as given in (^). The appearance of the logarithm of the determinant 

A = cxp Q Tr In [2m Q +iv-D-irj\j (48) 

is a quantum effect not present in the classical derivation presented above. 
However, in this case the determinant can be regulated in a gauge-invariant 
way, and by choosing the gauge v ■ A — one can show that In A is just an 
irrelevant constant EjEiJ. 
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Because of the phase factor in (fff|), the x dependence of the effective 
heavy-quark field h v is weak. In momentum space, derivatives acting on h v 
correspond to powers of the residual momentum k 7 which by construction is 
much smaller than vtiq. Hence, the non-local effective Lagrangian ( |46|) allows 
for a derivative expansion in powers of iD/niQ: 



£ e ff = h v iv-Dh v 



oo 

— y k 



iv ■ D 



2m 



Q 



ip ± h v . (49) 



Taking into account that h v contains a -P+ projection operator, and using the 
identity 



P+ip ± ip ± P+=P+ 



(iD±) 2 + f <V G^ 



+ ■ 



(50) 



where [i-D M , iD v ] — ig s G^ v is the gluon field-strength tensor, one finds thatBEl 

£ cff = h v iv-Dh v + -—h v (iD x ) 2 h v + -^-h v <Tav G" v h v + 0{l/m%) . (51) 
2m Q Aitiq * 

In the limit mq — > oo, only the first terms remains: 

Coo = h v iv-D h v . (52) 

This is the effective Lagrangian of the HQET. It gives rise to the Feynman 
rules depicted in Fig. ^. 



• k + in 2 J* 



ig S v a (t a )ji 



Figure 7: Feynman rules of the HQET and a are colour indices). A heavy quark is 
represented by a double line labelled by the velocity v and the residual momentum k. The 
velocity v is conserved by the strong interactions. 

Let us take a moment to study the symmetries of this LagrangianS Since 
there appear no Dirac matrices, interactions of the heavy quark with gluons 
leave its spin unchanged. Associated with this is an SU(2) symmetry group, 
under which Coo is invariant. The action of this symmetry on the heavy-quark 
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fields becomes most transparent in the rest frame, where the generators S l of 
SU(2) can be chosen as 



1 a 1 

2 V a 1 



[S i ,S j ]=ie ijk S k . (53) 



Here <j % are the Pauli matrices. An infinitesimal SU(2) transformation h v — > 
(1 + it - S) h v leaves the Lagrangian invariant: 

5£oo — h v [iv-D, ie ■ S] h v = . (54) 

Another symmetry of the HQET arises since the mass of the heavy quark does 
not appear in the effective Lagrangian. For Nh heavy quarks moving at the 
same velocity, eq. ( [52] ) can be extended by writing 

N h 

C 00 =^K i v ivDh i v . (55) 

i=l 

This is invariant under rotations in flavour space. When combined with the 
spin symmetry, the symmetry groupuis. promoted to SU(27V^). This is the 
heavy-quark spin-flavour symmetry F U 2 i. Its physical content is that, in the 
limit m,Q — > oo, the strong interactions of a heavy quark become independent 
of its mass and spin. 

Consider now the operators appearing at order 1/mg in the effective La- 
grangian (|5l|). They are easiest to identify in the rest frame. The first operator, 

Okin = K (*D±) 2 hy -> hy {lD f hy , (56) 

2mq 2niQ 

is the gauge-covariant extension of the kinetic energy arising from the off- 
shell residual motion of the heavy quark. The second operator is the non- 
abelian analogue of the Pauli interaction, which describes the chromo-magnetic 
coupling of the heavy-quark spin to the gluon field: 

Cmag — "; h v (Jul/ h v ► h v S • B c h v . (57) 

4mQ tuq 

Here S is the spin operator defined in and B\ — — jt l ^ k G^ k are the 

components of the chromo-magnetic field. The chromo-magnetic interaction 
is a relativistic effect, which scales like I/tuq. This is the origin of the heavy- 
quark spin symmetry. 
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3.2 The Residual Mass Term and the Definition of the Heavy-Quark Mass 

The choice of the expansion parameter in the HQET, i.e. the definition of the 
heavy-quark mass mQ, deserves some comments. In the derivation presented 
earlier in this section, we chose toq to be the "mass in the Lagrangian" , and 
using this parameter in the phase redefinition in ( [4l| ) we obtained the effective 
Lagrangian (|52|), in which the heavy-quark mass no longer appears. However, 
this treatment has its subtleties. The symmetries of the HQET allow a "resid- 
ual mass" 5m for the heavy quark, provided that 5m is of order Aqcd and is 
the same for all heavy-quark flavours. Even if we arrange that such a mass 
term is not present at the tree level, it will in general be induced by quantum 
corrections. (This is unavoidable if the theory is regulated with a dimensionful 
cutoff.) Therefore, instead of ( |5^ ) we should write the effective Lagrangian in 
the more general formEi: 

h v (x) = e im * v - x P + Q(x) 
=>■ Coo — h v iv-Dh v — 5mh v h v . (58) 

If we redefine the expansion parameter according to mQ — > tuq + Am, the 
residual mass changes in the opposite way: 5m — > 5m — Am. This implies that 
there is a unique choice of the expansion parameter mQ such that 5m = 0. 
Requiring 5m = 0, as it is usually done implicitly in the HQET, defineSpa-heavy- 
quark mass, which in perturbation theory coincides with the pole massEj. This, 
in turn, defines for each heavy hadron a parameter A (sometimes called the 
"binding energy" ) through 

A = (m H - m Q ) . (59) 

m,Q — >oo 

If one prefers to work with another choice of the expansion parameter, the 
values of non-perturbative parameters such as A change, but at the same time 
one has to include the residual mass term in the HQET Lagrangian. It can be 
shown that the various parameters that depend on the definition of mQ enter 
the predictions for all physical observables in such_a way that the results are 
independent of which particular choice one adopts l3. 

There is one more subtlety hidden in the above discussion. The quantities 
mQ, A and 5m are non-perturbative parameters of the HQET, which have a 
similar status as the vacuum condensates in QCD phenomenology Ej. These 
parameters cannot be defined unambiguously in perturbation theory. The rea- 
son lies in the divergent behaviour of perturbative expansions in large orders, 
which is associated with the existence of-ainftularities along the real axis in the 



Borel plane, the so-called renormalons E3 c3. For instance, the perturbation 
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series which relates the pole mass rriQ of a heavy quark to its bare mass, 
m Q = TO Q 1C I 1 + c i a s{m Q ) + c 2 a 2 s {m Q ) + . . . + c n a"(m Q ) + . .. j , (60) 

contains numerical coefficients c„ that grow .as n! for large n, rendering the 
series divergent and not Borel summableOO. The best one can achieve is to 
truncate the perturbation series at the minimal term, but this leads to an un- 
avoidable arbitrariness of order Atoq ~ Aqcd (the size of the minimal term) . 
This observation, which at first sight seems a serious problem for QCD phe- 
nomenology, should actually not come as a surprise. We know that because of 
confinement quarks do not appear as physical states in nature. Hence, there is 
no unique way to define their on-shell properties such as a pole mass. In view of 
this, it is actually remarkable that QCD perturbation theory "knows" about its 
incompleteness and indicates, through the appearance of renormalon singular- 
ities, the presence of non-perturbative effects. We must first specify a scheme 
how to truncate the QCD perturbation series before non-perturbative state- 
ments such as 8m = become meaningful, and hence before non-perturbative 
parameters such as tuq and A become well-defined quantities. The actual 
values of these parameters will depend on this scheme. 

We stress that the "renormalon ambiguities" are not a conceptual problem 
for the heavy-quark expansion. In fact, it can be shown quite generally that 
these ambiguities cancel in all predictions for physical observablest3~Ell. The 
way the cancellations occur is intricate, however. The generic structure of the 
heavy-quark expansion for an observable is of the form: 

observable ~ C[a s (m Q )} [ 1 + — + . • • | ■ (61) 

V m Q J 

Here C[a s (mg)] represents a perturbative coefficient function, and A is a 
dimensionful non-perturbative parameter. The truncation of the perturba- 
tion series defining the coefficient function leads to an arbitrariness of order 
Aqcd/toq, which precisely cancels against a corresponding arbitrariness of 
order Aqcd in the definition of the non-perturbative parameter A. 

The renormalon problem poses itself when one imagines to apply pertur- 
bation theory in very high orders. In practise, the perturbative coefficients are 
known to finite order in a s (at best to two-loop accuracy) , and to be consistent 
one should use them in connection with the pole mass (and A etc.) defined to 
the same order. For completeness, we note that for the "one-loop pole masses" 
of the heavy quarks we shall adopt the values 

m b = 4.8 GeV , m c = 1.4 GeV . (62) 

Their difference satisfies the constraint in (|l5|). 
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4 Inclusive Decay Rates 

Inclusive decay rates determine the probability of the decay of a particle into 
the sum of all possible final states with a given set of quantum numbers. An 
example is provided by the inclusive semileptonic decay rate of the B meson, 
T(B — > X c lv), where the final state consists of a lepton-neutrino pair accom- 
panied by any number of hadrons with total charm-quark number n c = 1. Here 
we shall discuss the theoretical description of inclusive decays of hadrons con- 



taining a heavy quark l§HI23. From the theoretical point of view, such decays 
have two advantages: first, bound-state effects related to the initial state (such 
as the "Fermi motion" of the heavy quark inside the hadron) can be accounted 
for in a systematic way using the heavy-quark expansion, in much the same 
way as explained in the previous sections; secondly, the fact that the final state 
consists of a sum over many hadronic channels eliminates bound-state effects 
related to the properties of individual hadrons. This second feature is based on 
a hypothesis known as quark-hadron duality, which is an important concept 
in QCD phenomenology. The assumption of duality is that cross sections and 
decay rates, which are defined in the physical region (i.e. the region of time-like 
momenta), are cakmlable in QCD after a "smearing" or "averaging" procedure 
has been applied E23. In semileptonic decays, it is the integration over the 
lepton and neutrino phase space that provides a "smearing" over the invariant 
hadronic mass of the final state (so-called "global" duality). For non-leptonic 
decays, on the other hand, the total hadronic mass is fixed, and it is only the 
fact that one sums over many hadronic states that provides an "averaging" 
(so-called "local" duality) . Clearly, local duality is a stronger assumption than 
global duality. It is important to stress that quark-hadron duality cannot yet 
be derived from first principles, although it is a necessary assumption for many 
applications of QCD. The validity of global duality has been tested experimen- 
tally using data on hadronic r decays EH. A more formal attempt to address 
the problem of quark-hadron duality can be found in Ref. 108. 

Using the optical theorem, the inclusive decay width of a hadron H b con- 
taining a b quark can be written in the form 

T(H b ^X) = -^—21m(H b \T\H b ), (63) 
2m Hb 

where the transition operator T is given by 

T = ifd 4 xT{C eB (x),Ces(0)}- (64) 
In fact, inserting a complete set of states inside the time-ordered product, we 
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recover the standard expression 

T(H b ^X) = -^— y2(27r)H i (p H - Px )\(X\C cS \H b )\ 2 (65) 
2mH » x 

for the decay rate. HsceJ[Lffi is the effective weak Lagrangian corrected for 



short-distance effectsEZIO'EHHUlHI arising from the exchange of gluons with vir- 
tualities between mw and nib- If some quantum numbers of the final states 
X are specified, the sum over intermediate states is to be restricted appropri- 
ately. In the case of the inclusive semileptonic decay rate, for instance, the 
sum would include only those states X containing a lepton-neutrino pair. 
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Figure 8: Perturbative contributions to the transition operator T (left), and the correspond- 
ing operators in the OPE (right). The open squares represent a four-fermion interaction of 
the effective Lagrangian C e g, while the black circles represent local operators in the OPE. 

In perturbation theory, some contributions to the transition operator are 
given by the two-loop diagrams shown on the left-hand side in Fig. || Because 
of the large mass of the b quark, the momenta flowing through the internal 
lines in these diagrams are large. It is thus possible to construct an OPE for 
the transition operator, in which T is represented as a series of local operators 
containing the heavy-quark fields. The operator with the lowest dimension, 
d = 3, is bb. It arises from integrating over the internal lines in the first diagram 
shown in the figure. The only gauge-invariant operator with dimension d = 4 
is bipb; however, the equation of motion implies that between physical states 
this operator can be replaced by rribbb. The first operator that is different from 
bb has dimension d — 5 and contains the gluon field. It is given by b gga^G^b. 
This operator arises from diagrams in which a gluon is emitted from one of the 
internal lines, such as the second diagram shown in the figure. 

For dimensional reasons, the matrix elements of such higher-dimensional 
operators are suppressed by inverse powers of the heavy-quark mass. Thus, 
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any inclusive decay rate of a hadron H b can be written in the form! 



T(H b ^X f ) = 1 ^^4(bb) H + ci + ... j, (66) 

where the prefactor arises naturally from the loop integrations, c£ are calcu- 
lable coefficient functions (which also contain the relevant CKM matrix ele- 
ments) depending on the quantum numbers / of the final state, and (0)h arc 
the (normalized) forward matrix elements of local operators, for which we use 
the short-hand notation 

(0) H = 7 ^—(H b \0\H b ). (67) 

In the next step, these matrix elements are systematically expanded in 
powers of l/m&, using the technology of the HQET. Introducing the velocity- 
dependent fiplHft of the HQET, where v denotes the velocity of the hadron 
H b , one finds 



1 fil(H b )-^ G (H b 



(bb) H = 1 - ™"°> + 0(l/mi) , 

(bg s a^G^b) H = 2n 2 G (H b ) + 0(l/m b ) , (68) 
where we have defined the HQET matrix elements 

A(H b ) = (H b (v)\ k, {iD) 2 b v \H b (v)} , 

uliHb) = tt^— (H b (v)\ h^-a^G^K \H b (v)) . (69) 
Zm Hb I 

Here {iD) 2 = (iv ■ D) 2 — (iD) 2 : in the rest frame, this is the square of the 
operator for the spatial momentum of the heavy quark. Inserting these results 
into (fitf), we obtain 



b J 



It is instructive to understand the appearance of the "kinetic energy" contri- 
bution fx"^, which is the gauge-covariant extension of the square of the 6-quark 
momentum inside the heavy hadron. This contribution is the field-theory ana- 
logue of the Lorentz factor (1 — v 2 ) 1 ^ 2 — 1 — p 2 /2m 2 , in accordance with the 
fact that the lifetime, r = for a moving particle increases due to time 
dilation. 
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The main result of the heavy-quark expansion for inclusive decay rates is 
that the free quark decay (i.e. the parton model) provides the first term in a 
systematic l/mj expansion, i.e. 

T(H b - X t ) = ^ 4 {l + 0{l/ml)} . (71) 

For dimensional reasons, the free-quark decay rate is proportional to the fifth 
power of the 6-quark mass. The non-perturbative corrections to this picture, 
which arise from bound-state effects inside the hadron Hi,, are suppressed 
by (at least) two powers of the heavy-quark mass, i.e. they are of relative 
order (Aqcd/™&) 2 - Note that the absence of first-order power corrections is a 
simple consequence of the equation of motion, as there is no independent gauge- 
invariant operator of dimension d = 4 that could appear in the OPE. The fact 
that bound-state effects in inclusive decays are strongly suppressed explains a 
posteriori the success of the parton model in describing such processes. 

The hadronic matrix elements appearing in the heavy-quark expansion 
( fzoj ) can be determined to some extent fmni the known masses of heavy hadron 
states. For the B meson, one finds thatu 

H%{B) = -Ai = (0.3 ± 0.2) GeV 2 , 

fi%{B) = 3A 2 = | (m|. - m|) ~ 0.36 GeV 2 , (72) 

where Ai and A2 are the parameters appearing in the mass formula (|J), and 
their numerical values have been taken from (|l4|) and (||). For the ground-state 
baryon A&, in which the light constituents have total spin zero, it follows that 

Mg(A 6 ) = 0, (73) 

while the matrix element /u£.(A(,) obeys the relation 

{m Ab -m Ac ) - (m B -fn D ) = /i 2 (B) - & (A 6 ) - ^-^j + 0{l/m 2 Q ) , 

(74) 

where ms and mo denote the spin-averaged masses introduced in connection 
with (p"3|). With the value of mA b given in (pr|), this leads to 

lil (B) - nl(K b ) = (0.01 ± 0.03) GeV 2 . (75) 

What remains to be calculated, then, is the coefficient functions c£ for a given 
inclusive decay channel. We shall now discuss two important applications of 
this general formalism. 
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4-1 Determination of \V c b\ from Inclusive Semileptonic Decays 

The extraction of \V c b\ from thepinclusivc semileptonic decay rate of the B 
meson is based on the expression E3H121I 

(76) 

where m& and m c are the poles mass of the b and c quarks (defined to a given 
order in perturbation theory), and f{x) and g(x) are phase-space functions: 

f{x) = l-8x 2 + 8x 6 -x 8 -12x 4 lnx 2 , (77) 




and g(x) is given elsewhere 1113. The theoretical uncertainties in this determi- 
nation of \V c b\ are quite different from those entering the analysis of exclusive 
decays. In particular, in inclusive decays there appear the quark masses rather 
than the meson masses. Moreover, the theoretical description relies on the as- 
sumption of global quark-hadron duality, which is not necessary for exclusive 
decays. 

A careful analysis shows that the main sources of theoretical uncertainties 
are the dependence on the heavy-quark masses and unknown higher-order per- 
turbative corrections E13. Each lead to an uncertainty of 5T/T ~ 10% in the 
prediction for the semileptonic decay rate. Adding, as previously, the theoret- 
ical errors linearly, and taking the square root, leads to 

S\V cb \ 



\Vcb\ 



10% (78) 



for the theoretical uncertainty in the determination of \V c b\ from inclusive 
decays, keeping in mind that this method relies in addition on the assumption 
of global quark-hadron duality. Taking the result of Ball et al. for the central 
value E3, we quote 

-DSL \ I T B 



\V cb \ = (0.0400 ± 0.0040) (79) 

With the new world averages for the semileptonic branching ratio_.i?sL = 
(10.90 ± 0.46)% (see below), and for the average i?-meson lifetime Lj, tb = 
(1.60 ± 0.03) ps, we obtain 

\V cb \ = (40.0 ± 0.9 cxp ± 4.0 th ) x 10- 3 . (80) 
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This is in excellent agreement with the value in ( |37| ) , which has been extracted 
from the analysis of the exclusive decay B — > D*l v. This agreement is grat- 
ifying given the differences of the methods used, and it provides an indirect 
test of global quark-hadron duality. Combining the two measurements gives 
the final result 

\V cb \ = 0.039 ±0.002. (81) 
After V u d and V us , this is now the third-best known entry in the CKM matrix. 



4-2 Semileptonic Branching Ratio and Charm Counting 
The semileptonic branching ratio of B mesons is defined as 

d T(B-+XeD) 

SL E*r(s^x^) + r NL + r rare ' 1 ' 

where Tnl and r raro arc the inclusive rates for non-lcptonic and rare decays, 
respectively. The main difficulty in calculating Bql is not in the semileptonic 
width, but in the non-leptonic one. As mentioned above, the calculation of non- 
lcptonic decays in the heavy-quark expansion relies on the strong assumption 
of local quark-hadron duality. 

Measurements of the semileptonic branching ratio have been performed by 
various experimental groups, using both model-dependent and model-indepen- 



dent analysestLM The status of the results is controversial, as there is a discrep- 
ancy between low-energy measurements performed at the T(4s) resonance and 
high-energy measurements performed at the Z° resonance. The average value 



at low energies isEL£l £?sl = (10.37 ± 0.30)%. High-energy measurements per- 



formed at LEP, on the other hand, givelLia B^ = (11.11 ± 0.23)%. The super- 
script (b) indicates that this value refers not to the B meson, but to a mixture 
of b hadrons (approximately 40% B~, 40% B°, 12% B Sl and 8% A 6 ). Assuming 
that the corresponding semileptonic width is close to that of the B meson,[] 

we can correct for this and find B SL = (r(S)/r(&)) = (11.30 ± 0.26)%, 
where r(b) = (1.57±0J(13) ps is the average lifetime corresponding to the above 
mixture of b hadrons c3. The discrepancy between the low- and high-energy 
measurements of the semileptonic branching ratio is therefore larger than three 
standard deviations. If we take the average and inflate the error to account 
for this fact, we obtain 

B SL = (10.90 ±0.46)%. (83) 



Theoretically, this is expected to be a very good approximation. 
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In understanding this result, an important aspect is charm counting, i.e. the 
measurement of the average number n c of charm hadrons produced per B 
decay. Theoretically, this quantity is given by 



1 + B(B -> X c 



B(B — > no charm) , 



(84) 



where B(B — > X cSs i) is the branching ratio for decays into final states contain- 
ing two charm quarks, and B(B — ► po charm) ~ 0.02 is the Standard Model 



branching ratio for charmless decaysE£D~ELj. Recently, two new measurements 
of the average charm c ontent have been performed. The CLEO Collaboration 
ie@@ n c = 1.16±0.05, and the ALEPH Collaboration 
1.23 ±0.07. The average is 



has presented the value, 
has reported the result 
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= 1.18 ± 0.04. 



(85) 



The naive parton model predicts that i?sL — 15% and n c ~ 1.2; however, it 
has been known for some-time that perturbative corrections could change these 
predictions significantly till. With the establishment of the 1/mQ expansion, 
the non-perturbative corrections to the parton model could be computed, and 
their effect turned out to be very small. This led Bigi et al. to conclude that 
values £?sl < 12.5% cannot be accommodated by theory, thus giving rise to 
a puzzle referred to as the "baffling semileptonic branching ratio" t23. Later, 
Bagan et al. have completed the calculation of the 0(a s ) corrections including 



the effects o: 
significantly 



the charm-quark mass, finding that they lower the value of -Bsl 



1 2:; 



Ehe original analysis of Bagan et al. has recently been corrected in an erra- 



tumE£3. Here we shall present the results of an independent numerical analysis 
using the same theoretical input (for a detailed discussion, see Ref. 124). The 
semileptonic branching ratio and n c depend on the quark-mass ratio m c /mb 
and on the ratio /Lt/m&, where \i is the scale used to renormalize the coupling 
constant a s (/i) and the Wilson coefficients appearing in the non-leptonic decay 
rate. The freedom in choosing the scale [i reflects our ignorance of higher-order 
corrections, which are neglected when the perturbative expansion is truncated 
at order a s . Below we shall consider several choices for the renormalization 
scale. We allow the pole masses of the heavy quarks to vary in the range [see 
(§) and (H)] 



m b = (4.8 ±0.2) GeV, 



nib 



= (3.40 ±0.06) GeV, 



(86) 



corresponding to 0.25 < m c /m(, < 0.33. Non-perturbative effects appearing at 
order 1 jm^ in the heavy-quark expansion are described by the single parameter 



A2 — 0.12 GeV defined in (|72[); the dependence on the parameter Ai is the 
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same for all inclusive decay rates and cancels out in Bsl and n c . For the two 
choices fi = nib and fi = rrib/2, we obtain 

f 12.0 ± 1.0%; /x = "16, 
SL _ \ 10.9 ± 1.0%; /x = TO fc /2, 

_ f 1.20 T 0.06; /j = m b , 
Uc \ 1.21 T 0.06; fi = m b /2. (b '> 

The uncertainties in the two quantities, which result from the variation of 
m c / m b in the range given above, are anticorrelated. Notice that the semilep- 
tonic branching ratio has a stronger scale dependence than n c . By choosing 
a low renormalization scale, values Bsl < 12% can easily be accommodated. 
The experimental data prefer a scale /j/mui' 0.5, which is indeed not unnat- 
ural. Using the BLM scale setting method E3, \"\ p et al. have estimated that 
/i > 0. 32to;, is an appropriate scale in this caseE£3. 

1.4 



1.3 



c° 1.2 



1.1 



1 

8 9 10 11 12 13 14 
B SL (%) 

Figure 9: Combined theoretical predictions for the semileptonic branching ratio and charm 
counting as a function of the quark-mass ratio m c /mi, and the renormalization scale fi. The 
data points show the average experimental values for Bsl and n c obtained in low-energy 
(LE) and high-energy (HE) measurements, as discussed in the text. 

The combined theoretical predictions for the semileptonic branching ratio 
and charm counting are shown in Fig. H. They are compared with the ex- 
perimental results obtained from low- and high-energy measurements. It was 
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argued that the combination of a low semileptonic branching ratio and a law 
value of n c would constitute a potential problem for the Standard Model tlB. 
However, with the new experimental and theoretical numbers, only for the 
low-energy measurements a small discrepancy remains between theory and ex- 
periment. Note that, using (jSJ), our results for n c can be used to obtain a 
prediction for the branching ratio B(B — > X c£s >), which is accessible to a direct 
experimental determination. Our prediction of (22 ± 6)% for this branching 
ratio agrees well with the preliminary result reported by the CLEO Collabo- 



rationEfl: B(B -> X c5s ,) = (23.9 ± 3.8)%. 



5 Concluding Remarks 

We have presented an introduction to the theory and phenomenology of heavy- 
flavour physics. The theoretical tools that allow us to perform quantitative 
calculations in this area are heavy-quark symmetry, the heavy-quark effective 
theory, and the 1/m.Q expansion. After presenting the underlying theoretical 
concepts, we have discussed their application to the description of exclusive and 
inclusive weak decays of B mesons. Besides presenting the status of the latest 
developments, our hope was to convince the reader that heavy-flavour physics is 
a rich and diverse area of research, which is at present characterized by a fruitful 
interplay between theory and experiments. This has led to many significant 
discoveries and developments on both sides. Heavy-quark physics has the 
potential to determine many important parameters of the electroweak theory 
and to test the Standard Model at low energies. At the same time, it provides 
an ideal laboratory to study the nature of non-perturbative phenomena in 
QCD, still one of the least understood properties of the Standard Model. 

Let us finish with a somewhat philosophical remark: In the last few years, 
there have been exciting developments in high-energy physics related to the dis- 
covery dualities, which relate apparently very different theories to each other. 
So are electric, weak-coupling phenomena in one theory dual to magnetic, 
strong-coupling phenomena in another theory. Some people argue quite con- 
vincingly that duality seems to be everywhere in nature, and consequently 
there are no really difficult questions in physics; very difficult problems be- 
come trivial when approached from a different, dual point of view. There are, 
however, "moderately difficult" problems in physics, which are "self-dual" . 
It is the author's opinion that real-world (i.e. non-supersymmetric) QCD at 
hadronic energies belongs to this category. Having said this, we conclude that 
heavy-quark effective theory provides a powerful tool to tackle the "moderately 
difficult" problems of heavy-flavour physics. 
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